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We study the effects of hole doping on the two-dimensional Heisenberg-Kondo model around 
the quantum critical point, where the spin liquid phase (Kondo insulator) and the magnetically 
ordered phase are separated via a second-order phase transition. By means of the self-consistent 
Born approximation within the bond operator formalism as well as the standard spin wave 
theory, we discuss dynamical properties of a doped hole. It is clarified that a quasi-particle 
state stabilized in the spin liquid phase is gradually obscured as the system approaches the 
quantum critical point. This is also the case for the magnetically ordered phase. We argue the 
similarity and the difference between these two cases. 
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1. Introduction 

Heavy fermion systems realized in rare-earth com- 
pounds have attracted much interest over many years. 
They show a variety of striking phenomena due to strong 
electron correlations. In metallic systems, the electron 
mass is renormalized to a huge value by the /-/ Coulomb 
interaction. 1 - 1 Such renormalization effects are also rele- 
vant for the insulating phase, leading to the Kondo insu- 
lator possessing both of the spin and charge gaps renor- 
malized by electron correlations. 2 ' 3 ' Since these metal- 
lic and insulating states have essentially the same origin 
of electron correlations, it is interesting to ask how the 
Kondo insulator is changed to the heavy fermion metallic 
system upon hole doping. 

Another interesting aspect in heavy fermion systems 
is that some of unusual static and dynamical proper- 
ties are related to the fact that the system is located 
in the vicinity of a quantum critical point. 4 ~ 14 ' For in- 
stance, the anomalous non-Fermi-liquid temperature de- 
pendence in the specific heat and the susceptibility ob- 
served in rare-earth compounds such as CeCu6- x R x (R = 
Au, Ag), 15) CePd 2 Si 2 , 16 ) CeNi 2 Ge 2 , 17 ) Ui_ x Y x Pd 3 , 18 ) 
Ce x Lai_ x Ru 2 Si 2 19 - ) is caused by large quantum fluctua- 
tions reflecting the phase transition to the antiferromag- 
netically ordered phase. In this way, strong quantum 
fluctuations near the critical point 20 ) give rise to rich 
and attractive phenomena in rare-earth systems. 

Stimulated by these interesting topics, we consider a 
simplified model of the Kondo insulator, which possesses 
the quantum critical point between the Kondo insulat- 
ing phase and the antiferromagnetic phase. To be pre- 
cise, we employ the two-dimensional Heisenberg-Kondo 
model, for which the charge degree of freedom is frozen at 
half filling, and the remaining spin sector is either in the 
spin liquid phase or the antiferromagnetically ordered 
phase at zero temperature. We investigate the effects of 
hole doping on the model with particular emphasis on the 
hole dynamics near the quantum critical point. Although 
the effects of hole doping into the Kondo insulator have 
been addressed so far, 3 ) we think that the dynamics of 
a doped hole has not been discussed in detail. Also, in 
some respects, the present investigation is closely related 



to those done for the two-dimensional Heisenberg model 
in connection with the high T c superconductors. 21 ~ 24 ) 

This paper is organized as follows. In Sec. 2, we intro- 
duce the two-dimensional Heisenberg-Kondo model and 
summarize the basic properties of the model. In Sec. 3, 
we study the hole dynamics, when it is doped into the 
spin liquid phase, by combining the self-consistent Born 
approximation with the bond-operator formalism. Simi- 
lar analysis is performed for the antiferromagnetic phase 
in Sec. 4 by exploiting the spin wave theory. A brief 
summary is given in Sec. 5. 

2. Heisenberg-Kondo model 

We investigate the two-dimensional Heisenberg-Kondo 
model having conduction electrons coupled to localized 
/-spins. The model is schematically drawn in Fig. 1. 
The Hamiltonian reads 




Fig. 1. Heisenberg-Kondo model on the square lattice. Spins of 
host electrons indicated by the arrows are antiferromagnetically 
coupled to localized /-spins indicated by the circles. 



H = H t +Hj, (1) 
Hj = J 1 1 S?.S° + J ± £S«.S*, (2) 



<hj> 



H t = -t 4,a £ j,a+H.C. 



(3) 



<i,j>,a 



with Ci tC! = [cj )Cr (l — n^-o-)], where Cj iCr annihilates a 
conduction electron with spin a{=],{) at the ith site. 
The corresponding spin operator is defined by = 
1Eoj' c !j t m'Vi where r is the Pauli matrix. The 
operator of an s = 1/2 /-spin is denoted by Sj. The 
exchange couplings J// and Jj_ are assumed to be an- 
tiferromagnetic. For convenience, we set hopping t as 
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unity and introduce the ratio j = J///J±. At half fill- 
ing, the system is in the spin-liquid insulating phase with 
both of the spin and charge gaps for smaller j, which may 
be a proper model of the Kondo insulator. We note that 
magnetic properties for j = have already been dis- 
cussed. 25 ~ 29 ) An interesting aspect in the model is that 
there exists a quantum phase transition at j c (w 0.75, 8- 
th order series expansion) 30 ) from the spin-liquid phase 
(Kondo insulator) to the magnetically ordered phase as 
j increases. Therefore, the system realizes a quantum 
critical point at j — j c , where the spin gap vanishes, but 
a long-range order still does not develop. This specific 
condition may give rise to anomalous behavior of a doped 
hole, reflecting large quantum fluctuations. 

In order to discuss dynamical properties of a hole- 
doped Kondo insulator in the vicinity of the quantum 
critical point, we make use of either the bond-operator 
approach 31-33 ) or the spin wave approach. These meth- 
ods are appropriate to discuss magnetic properties in the 
spin liquid phase and the magnetically ordered phase, 
respectively. We further exploit the self-consistent Born 
approximation(SCBA), 21 ~ 24,34 ) in which the motion of 
hole can be treated as a scattering problem in the host 
spin system. In the following, we separately discuss dy- 
namical properties of a doped hole in the spin liquid 
phase and the magnetically ordered phase. 

3. Hole doping in the spin-liquid phase 

We begin with dynamical properties of a doped hole in 
the spin liquid phase (Kondo insulator) having the spin 
gap for excitations. For this purpose, we use the bond- 
operator representation, which is useful to deal with the 
spin-singlet state. 

3.1 bond-operator representation 

When the Kondo coupling between the host spins and 
localized spins is fairly strong, the dimer singlet is formed 
at each site, which may be a good starting point to de- 
scribe the spin liquid phase. This condition allows us 
to introduce six kinds of the bond operators, which are 
defined at each site as, 



4lo> 
4„l<» = 



IT)), 

m-iu)), 



4,jo) = ^=(im + iu», 

4,J0) - ^(ITD + IIT)), 



Vt 



Vi 



10) = 
10) = 



l°T>, 
1° i>, 



(4) 



where ket states are specified by the configuration of 
a conduction electron and an /-electron [o represents 
an unoccupied (hole) site]. Here, the operators st and 
ta( a = x tVt z ) obey the bosonic commutation relation, 
while the operator the fermionic anticommutation re- 



lation. To restrict the physical states to either singlet, 
triplet or pscudofcrmion, the above operators are sub- 
jected to the constraint, 



(5) 



In the spin liquid phase, the boson s is condensed, 
which results in the finite value of < s n >= s. Then the 
spin Hamiltonian H j is diagonalized in the Fourier space 
as, 



Hj = ^w k 7i ik 7a,k + const., 



(6) 



wk = JjVI + 2e k , 

where e k = js 2 Q(k) with Q(k) = ^(cosk x + cosfc y ). 
Here, we have introduced the normal-mode operator 7^ k 
related to t' a k via the Bogoliubov transformation, 



4k = u k7«,k + Vk7a,-k, 



(7) 



where w 2 = 1/2{J±(1 + e k )/co k - 1}, w 2 = 1/2{J±(1 + 
e k )/^ k + 1}. 

When a hole is doped into the system, the condensate 
density of the singlet is determined by the hard-core con- 
straint eq. (5) as 



(8) 



In the bond-operator representation, the operator for a 
physical electron is written by two bond operators as, 

+ 4, CT (P^x,n + ity,n)}, (9) 

where p a = +(— ), o =1 (t) for a =f (J.). We thus obtain 
the Hamiltonian for electron hopping, 

Ht = \f ]T a[ a a ha -ts ]T (ttTij+ttTj-O 

<i,j>,a <i,j> 

+1 Yl titA *,' - 1 Y. *(*< x ^) T ^ 

<ij> 

(10) 

with tj - (t J x<i ,tl.,tl.) and 
T = - V at 



<i,j>.a 

+H.c. 
i < - r+t ft f t 



2 / j -m,a 1 T o'i,o'2 a n,cr2j 



(11) 



where (T TOi „)t = T„, m . This Hamiltonian is spin rota- 
tionally invariant, reflecting the fact that the system has 
no magnetic order. 

3. 2 Green function 

To discuss the dynamics of a doped hole, we define 
the retarded Green functions for a pseudofermion and a 
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physical electron as 

G„(k,i) = -iQ(t)<D\{a Ka {t),aL ia }\D>,(12) 

G£(k,t) = - l e(t)<D\{c Ka (t),ctj\D>,(13) 

where \D > represents the spin-singlet ground state. For 
small j, rcnormalization effects due to the two triplet 
vertices in eq. (10) are of minor importance, so that we 
will discard them in the following. We then arrive at the 
Hamiltonian for a doped hole as 

H = 2ts 2 ^Q(k)a k(T a k , (T , (14) 

k 

(15) 



V 



k,q,(7i,(72 

+M kq 7_ q a[_ T ai :<J2 a k , CT2 } , 



where the element M in the perturbation term V repre- 
sents the scattering of a triplet boson, 



M kq = ts{Q(k) + Q(k - q)}(u q + w q ). 



(16) 



We wish to note that the SCBA is useful to deal with 
the Hamiltonian eqs. (14) and (15). In terms of the bare 
Green functions for a pseudofcrmion and a triplet-boson, 



G°(k,c) 



1 



u - 2ts 2 Q(k) ' 



ui — w k 



(17) 



(18) 



we can evaluate the self-energy of a pseudofcrmion in the 
SCBA, which is diagramatically shown in Fig. 2. The 

D (q ,CO q ) 



M. / 

K q i 



\M 



kq 



G a {k- q,(0 - (O q ) 

Fig. 2. Self-energy of a pseudofcrmion in the SCBA. The triangle 
represents the scattering matrix of a triplet boson. 

Green function defined by G CT (k, uj) = [G°(k, w)^ 1 — 
£<j(k, uj)]^ 1 is determined via the self-consistent equa- 
tion, 

S CT (k,c) = l^TM^G^k-q^-^)- 1 



-S CT (k- q,w - uq)] \ 



(19) 



where u) = w + ir/ with ry — > + . 

To obtain the physical Green function of an electron, 
we introduce the elements as 

±s 

ak, CT = s(D\a Ks c Ka \D) = 

/3kq,o- = (£>|(7q a k-q)(l/2)(Tl/2)Ck,cr|£>) 

T 



(20) 



which relate the Green function of a pseudofermion to 
that of a real electron. When j < j c , the density of 
triplets is small in the ground state, i.e., < t), < 1. 
Then it is legitimate to use a perturbative evaluation 
of two-particle propagator within the random-phase ap- 
proximation (RPA). 34 ) This gives us the Green function, 

G%{k,uj) = a k ^G rT (k,uj)a k ^ 

+ ^f3kq,aG a (k - q, OJ - UJ q )f3 k q.a 



+ (a k ^G a (k,uj),^^^2p kq ^G a (k -q,u- w 9 )M^ 



1 

-G a {k,w) 



%E g M kg G a (k-q,oj-oj g )Ml 



kg 



^E g M fcgG CT (fc - q,uj -UJ q )P kq ^ \ 
\ G rT (k,uj)a k , a J 

This completes the formulation of the SCBA for a hole 
doped in the spin liquid phase. 

3.3 numerical results 

Before discussing the nature of the hole-propagator, 
we check how well our bond-operator approach works at 
half filling. In Fig. 3, we show the spin-triplet excita- 
tion gap calculated by several different methods. The 
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Fig. 3. Spin-triplet excitation gap for the Heisenberg-Kondo 
model at half filing. The solid, dash-dotted, and broken lines 
represent the spin gap computed by the bond-operator method, 
the series expansion method (9-th order), and the exact diago- 
nalization (3x4 system). 

phase transition point was evaluated approximately as 
j c ~ 0.75, 30 ) so that we can say that the bond-operator 
approach may describe the spin-liquid phase rather well 
in a wide range of the parameters. 

Let us now turn to the hole-doped case. From now on, 
we set i] = 0.1. Shown in Fig. 4 is the spectral function 
of a doped hole computed by the SCBA. When J± = 5 
and j — 0.5, we find a dominant dispersive peak in the 
spectral function, which originates from the first term 
in the Hamiltonian (10), indicating the coherent motion 
of a doped hole without triplet excitations. Namely, the 
low-energy edge of the spectrum determined by the poles 
of the Green function on the real axis remains separated 
by gap from the continuum of multimagnon shake off. 
The corresponding dispersion of a quasiparticle is given 
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5 

(fl/t w/t 

Fig. 4. Spectral function of a doped hole in the spin liq- 
uid phase for J± = 5, j = 0.1 (a) and j = 0.5 (b). 
The momentum of each spectral function is chosen as k = 
(0, 0), (tt/4, tt/4), (it/2, tt/2), (3tt/4, 3tt/4), (it, it) from the bot- 
tom to the top. 



by 



e k = 2^Q ( k)--]T o 
q 



(22) 



Shown in Fig. 5 is the quasi-particle dispersion thus ob- 





j=o.i 







(-71, -71) 



(7U, 7U) 



Fig. 5. Dispersion relation of a doped hole for J± = 5 and j = 
0.1. The solid (broken) line is obtained by the SCBA (scries 
expansion of the 7th order). 

tained. The comparison of this with the spectral func- 
tion in Fig. 4 confirms that the sharp peak is indeed 
developed along the dispersion curve. This implies that 
the quasi-particle state is quite stable in this parame- 
ter region. We also calculate the dispersion of a hole by 
means of the series expansion, 35 ) in which the excitation 
energy is expanded in j and t up to the seventh order. 
The results are shown in Fig. 5 by the broken line. We 
find that the results obtained by the SCBA are consis- 
tent with those obtained by the series expansion. In- 
creasing the interdimer coupling J//, antiferromagnetic 
correlations are gradually enhanced and the spin gap is 
decreased. Then a doped hole suffers from the scattering 
by low-lying spin excitations, making the peak structure 
on the dispersion curve somehow broadened. This ten- 
dency is clearly observed in the peak around k = (0, 0) 
in Fig. 4(b). When J± = 1 and j = 0.5, the spin gap be- 
comes small. In this case, the quasiparticle peak around 
k = (0, 0) is completely smeared, as shown in Fig. 6. 
The large portion of the corresponding weight is shifted 
to the incoherent parts in higher energy region. The 
dispersion relation obtained approximately possesses a 
dip structure around the origin, as seen in Fig. 6(b). 
However, this structure may not be sensible, since the 
quasi-particle picture does not hold in this momentum 




(7T,7t) 



Fig. 6. Spectral function of a doped hole (a) and the corresponding 
dispersion relation (b) for J± = 1 and j = 0.5. The momenta 
are chosen as in Fig. 4. 



region. We note that even in this case, the spectrum 
around k = (it, n) still forms a rather sharp structure, 
implying that the hole motion with this momentum is 
not affected so much by spin excitations. For reference 
we show the width W of the quasi-particle band and the 
renormalization factor Z of the k = (0, 0) state in Fig. 
7. The latter quantity is defined by 



Z(k) 



<9£(k,w) 
duo 



(23) 



which corresponds to the weight of the quasi-particle 
state. It is seen that the band width becomes small as 
the system approaches the quantum critical point. At 
the same time, the weight of the quasi-particle decreases 
monotonically near the critical point, being consistent 
with the disappearance of the well-defined quasi-particle 
state. We have so far discussed spectral properties of a 
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Fig. 7. Plots of the width W of the quasi-particle band and the 
renormalization factor Z for the spin-liquid phase. 

doped hole in the spin liquid phase, and shown that well- 
defined quasi-particle states in the j <C j c are gradually 
obscured when the system approaches the quantum crit- 
ical point j = j c . In particular, the spectrum around 
k = (0, 0) is affected considerably. It is to be noted here 
that these characteristic properties may not be specific 
to the present model, but more generically hold for hole- 
doped systems in the spin liquid phase. 

Before closing the section, we wish to make closer com- 
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parison of the present SCBA with the exact diagonaliza- 
tion calculation by taking the one-dimensional model. 
In Fig. 8, we show the spectral function and the disper- 
sion relation for the one-dimensional Hciscnbcrg-Kondo 
model. It is found that the results of the SCBA agree 
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Fig. 8. (a): Spectral function of a doped hole in the onc- 
dimcnsional Kondo insulator at J±_ = 5. The solid line is ob- 
tained by the SCBA and the dashed line by the exact diago- 
nalization for the system of 10 sites. The momentum is varied 
as k = 0, 7r/4, tt/2, 37r/4, it. (b): The solid line is the dispersion 
relation obtained by the SCBA at J± = 5. The triangles and 
the squares represent the results obtained by the exact diago- 
nalization for the system of 8 and 10 sites, respectively. 

fairly well with those of the exact diagonalization, both of 
which show the well-defined quasi-particle state. In this 
way, even in one dimension, the existence of the spin gap 
plays an important role in stabilizing the quasi-particle 
state in the spin liquid phase. 

4. Hole doping in the antiferromagnetic phase 

We next discuss the effects of hole doping on the an- 
tiferromagnetically ordered phase j > j c . In this case, 
excitations in the host spin system may be well described 
by the spin wave theory, which we will use in the follow- 
ing discussions. 

4-1 spin-wave theory 

We exploit the Holstein-Primakoff transformation, 
which expresses the spin-1/2 operators Sf and in 
terms of the Bose operators 6 and d as, 



sr 



= b\J{l-b\b % ) 



b\b, 



4y/{l-4di) 
y/(l-dldi)di 
-\+4di 



(24) 



for each sublattice. By performing the Fourier trans- 
formation and then applying the standard Bogoliubov 
transformation, the Hamiltonian Hj for the spin sector 
can be diagonalized as 



Hj = ^(wika^ak + w 2k /3 k /3k) + const., 



(25) 



with cji k = 4J//VS k — A k and w 2 k = 4J//V_B k + A k , 
where A k = 1/16[4 + ij' 1 + Q(k) 2 (r 2 - 4.T 1 - 8) + 
4Q(k) 4 ] 1/2 , Bk = 1/8(1 - Q(k) 2 +j~ 1 /2). Here, we 
have ignored the hard-core constraint in cq. (24) since it 
has little effect on low-energy properties for small values 
of j -1 . The Bogoliubov transformation from the original 
bosons (6k) <^k) to the normal-mode bosons (a k , /3 k ) takes 
the form, 



6 k 



V dl k J 

where 



r k 
-e k 



-e k 
r k 



/ a k \ 
V ft* I 



(26) 



and 



r k 



8Sk, : 



^?k,+ 










Lk,- 


1 


1 



(27) 



r 



Pk,± 

2 



Q(k)F k ,± 
1 



{- 



64 



E^±(-+F K± )}/p K± 



Q(k)i" 



with p k ± = i? k ± + I£ ± - Li ± - 1. Here, E K ± = 

j-78 + VS k T Ak, F k> ± =j- 1 /8- V^kT A k . 

We show the dispersion relation obtained for the un- 
doped case in Fig. 9, in which gapless excitations in 





G) lk 






G) 2k 













(-7T,-7T) 



(71,71) 



Fig. 9. Spin wave spectra at J// = l,j = 0.25 along the line 
in the momentum space from k = (— tt, —it) to k = (7r,7r). The 
solid line indicates oj lk and the dashed line Wie- 
the low-energy region u>i and gapful excitations with a 
large dispersion w 2 appear. It is easily seen that u>i is 
mainly contributed by localized spins, while w 2 by con- 
duction electrons in this parameter region. To discuss 
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the stability of the magnetically ordered ground state, 
we calculate the spin deviation A Si for /-spins and AS2 
for conduction electrons, 



ASi =< 4dk >= ^£«+ +Pk,-), 

k 

A5 2 =< blb k > = 1 £(i£ + + L l-)- ( 29 ) 



In Fig. 10, we show the staggered magnetization S — AS 
as a function of j' -1 . When J± = = 0), the 




Fig. 10. Staggered magnetizations S — ASi (solid line) and S — 
AS2 (dashed line) as a function of 

system is reduced to the Heisenberg antiferromagnet on 
the square lattice, 36 - 1 which is completely separated from 
free /-spins S b . The introduction of the exchange cou- 
pling J± enhances spin correlations between conduction 
electrons and localized spins. Then the magnetization 
of conduction electrons once increases and has a max- 
imum value around j" 1 w 0.37. Further increase of 
the exchange coupling enhances the dimer correlation, 
and thereby suppresses the magnetization. Note that 
the ordered state is stabilized up to the critical value 
j^ 1 w 6.97 for conduction electrons, but j^ 1 w 7.60 for 
localized spins. This pathological result means that be- 
yond j~ x w 6.97 our spin-wave theory may break down. 
Since the critical value j~ x s=s 6.97 is slightly larger than 
that obtained by the series expansion, we believe that our 
spin wave theory can properly describe the ordered state 
except for the region very close to the phase transition 
point. 

4-2 Green function 

Let us now turn to the hole-doped case. In the spin- 
wave approximation, the Hamiltonian for electron hop- 
ping reads 



H, 



-* £ {^M^ 1 - 6 &) 



<*j> 



+ (1 - b\bj)bi] + H.c.}, 



(30) 



where /i,(= c\ a ) is the annihilation operator of a hole. 
This is rewritten in terms of the normal modes as, 



£ 4^k-q[(#q, + 7k-q - £q, + 7k)(a q + 
k,q 

+ (fl q ,-7k-q - iq,-7k)(/3q + Pla)]- 



We introduce the following Green function of a doped 
hole, 



G(M) = -iMT[Mt)4(o)M, 



(32) 



where \ip) is the ground state of the model and T is a 
time-ordering operator. When a doped hole hops in the 
lattice, the ordered state is perturbed and magnons are 
excited. Note that the motion of a hole scatters two 
types of magnons specified by c*k and /3k- By exploiting 
the SCBA for magnon scattering, we have the self-energy 
of the hole Green function as 



£(k,w) 



£«(k,w) + £( 2 )(k,w) 



(33) 



£«(k, W ) = ^^(i? q , +7 k-q-L q , +7k ) 2 G(k-q,c-c 1:q ) 



N 



£( 2 >(k,w) = ^£(i? q ^7k- q -i q ,-7k) 2 G(k-q,c-c 2 , q ). 



N 



When J ' 1 1 3> J_i_ and t, one obtains the dispersion 
relation of a hole analytically as, 



£k = 

£ik = 

6k = 



— £lk — £,2k 

m 2 ^ (i? q , + 7k-q - iq, + 7k) 2 



N 



£ 



Wlq 



16t 2 (flq,-7k-q - £q,-7k) 2 



N 



£ 



^2q 



• (35) 



^.3 numerical results 

We show the spectral function in Fig. 11(a). When a 
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Fig. 11. (a) Spectral function of a doped hole. The momenta are 
chosen as in Fig. 4. (b) Dispersion relation given by Eq. (35): 
j = 0.5, J ± = 10- 5 . 

hole is doped into the system, the peak structure devel- 
ops along the dispersion curve shown in Fig. 11(b), being 
consistent with those discussed in the doped Heisenberg 
antiferromagnet on the square lattice. 21 ~ 24 ) In partic- 
ular, the peak is rather sharp in the low-energy region 
near k w (ir/2, tt/2). These results imply that the quasi- 
particle state defined by eq. (35) is stable in this param- 
eter region. 

Shown in Fig. 12 is the spectral function A(k, u>) when 
the exchange coupling J± is varied. For j^ 1 = 0.25, the 
quasi-particle state is well defined, which has a visible 
•j dispersion in the spectral function. When the exchange 
~@oupling J± increases, the quasiparticle state becomes 
less dispersive. Such flattening effect is more clearly seen 
(33} direct plots of the dispersion relation shown in Fig. 



(34) 
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co/t co/t 



Fig. 12. Spectral function of a doped hole as a function of ui for 
J//=0.8. The momenta are chosen as in Fig. 4. 




k k 

Fig. 13. Dispersion relation (bold line) of a doped hole in the an- 
tiferromagnetic phase for J / / = 0.8. Thin and broken lines 
represent the contributions from £ lk and £ 2 k m eq.(35). 



13. The dispersive band for smaller J± is mainly con- 
tributed by conduction electrons. On the other hand, 
for larger J±, conduction electrons and /-electrons are 
mixed up, which in turn gives rise to the flat dispersion 
near the critical point. In this region, other peaks with 
large weight appear besides the quasiparticle peak in the 
spectral function, as seen in Fig. 12. This implies that 
the quasiparticle picture does not hold anymore. 

In Fig. 14 we show the width of the quasi-particle band 
and the renormalization factor. It is seen that both of 
these quantities become small as the system approaches 
the critical point, similarly to the case of the spin liq- 
uid phase discussed in the previous section. However, 
in closer comparison of these two cases, we see that the 
quasi-particle is less stable in the ordered phase, because 
the quasi-particle state in this case is dressed by gapless 
magnon excitations, making it somewhat difficult to dis- 
tinguish the quasi-particle state from low-energy excited 
continuum. 

5. Summary 

We have investigated dynamical properties of a sin- 
gle hole doped into the Hciscnbcrg-Kondo model on the 
square lattice. At half filling, the model is either in the 
disordered spin-liquid phase or the antiferromagnetically 
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Fig. 14. Band width W for the quasi-particle and renormalization 
factor Z in the ordered phase. We show Z with k = (it/2, it/2), 
since this state has a rather clear peak structure. Solid circles 
are the results for the 8x8 system. The solid line for the width 
is computed from the analytic expression eq.(35). It is not easy 
to extract the correct values of W from the spectral function 
computed for larger J±, so that we have plotted them only for 
the small region. The non-monotonic behavior in Z reflects 
a singular property of the limit J ± — > 0. 



ordered phase. In the spin liquid phase, by employing 
the bond-operator representation, we have shown that 
a doped hole has the nature of a well-defined quasi- 
particle with the unambiguous dispersion relation. In 
the vicinity of the quantum critical point, however, the 
quasi-particle state is obscured due to the decrease of 
the spin gap. This effect appears most prominently in 
the region around k = (0,0). Even in this case, it has 
shown that the spectrum around k = (jr, tt) still forms 
a sharp peak structure. In the antiferromagnetically or- 
dered phase, we have encountered similar quasi-particle 
behavior of a doped hole, as in the spin-liquid phase. 
Namely, the quasi-particle state is smeared and its dis- 
persion becomes flat as the system approaches the quan- 
tum critical point. 

Although apparent properties are quite analogous in 
these two cases, there is the essential difference between 
them for the origin of the quasi-particle state: the quasi- 
particle in the spin liquid phase is stabilized in the pres- 
ence of the spin gap, while that in the ordered phase is 
dressed by low-energy magnon excitations. We have in- 
deed seen that this causes the difference in the stability 
of the quasi-particle state around the critical point. 

In this paper, we have discussed the hole-doping effects 
by starting from two extreme limits. In these approaches, 
it is not easy to deal with the properties precisely just at 
the critical point. It is thus desirable to directly investi- 
gate anomalous properties such as the temperature de- 
pendence of the susceptibility around the critical regime. 
Also, it remains an interesting open problem what kind 
of unusual properties are expected when the finite den- 
sity of holes are doped into the system around the critical 
point, which is now under consideration. 
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